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Abstract 

The semileptonic and leptonic decays of heavy mesons are studied as a phenomeno- 
logical apphcation and exploration of a heavy-quark limit of Dyson-Schwinger equa- 
tions. The single form factor, ^(to), which characterises the semileptonic decay in 
this limit, is calculated and compares well with recent experimental extractions. 
We obtain a lower bound of 1/3 on the slope-parameter = —^'(1), which, in 
calculations with realistic input, is exceeded by a great deal: agreement with exper- 
imental data requiring ~ 1.2- 1.6. The flavour and momentum dependence of the 
light-quark propagators has observable consequences. 
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1. Introduction. 

Semileptonic decays of pseudoscalar mesons provide a means of measuring el- 
ements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix; fundamental pa- 
rameters of the standard model. For example: the K^.^ decay, K — >■ Trei^e, can 
be used to determine \Vus\] -Des to determine and semileptonic decays of 
B mesons can provide information about \Vcb\ and \ Vub\- With a single hadron 
in both the initial and final state, these processes are an ideal tool for studying 
the influence of nonperturbative, strong interaction dynamics on weak interac- 
tion processes. An accurate determination of the CKM matrix elements relies 
on developing a good understanding of these nonperturbative effects. 
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The simplest of these processes to study are those with a pseudoscalar meson 
in both the initial and final state. In this case the interaction is only sensitive 
to the vector part of the electroweak {V — A) interaction and, of the two form 
factors, denoted f±(t), f+(t) is dominant for e and /i final states, for which 
the leptonic current is approximately conserved. In theoretical studies of these 
decays of light pseudoscalar mesons, tt and K, nonperturbative dressing of the 
quark- W-boson vertex is important, with nonanalytic contributions associated 
with, for example, K-tt loops, being significant for t > [1]. For heavy mesons, 
e.g. D and B, such nonanalytic contributions are unimportant in the physical 
region: < t < (m^ — mo)^. However, this does not automatically mean that 
one can use the bare vertex, 7^, because these are not the only dressing effects 
induced by the strong interaction. 

The propagator for a dressed quark of fiavour / can be writtenQ 



The characteristic behaviour of Sf{p) in QCD can be estimated by studying 
the quark Dyson-Schwinger equation (DSE) [2] and is illustrated in Fig. 1. As 
typical of these calculations, the different quark flavours are specifled by the 
value of the current-quark mass at the renormalisation point, ~ 10 GeV in 
this heuristic illustration using the simple model of Ref. [3], with (in MeV): 
^uifi) ~ 1; ms{fi) ~ 25; rric^fi) ~ 300; and mfe(/x) ~ 3000. The intersection of 
the dotted, diagonal line with a given curve marks the solution of p = Mf{p) 
and deflnes the Euclidean constituent quark mass. The ratio Mf /mf 

is a single, indicative and quantitative measure of the nonperturbative effects 
of gluon-dressing on the quark propagator. 

Obvious in this flgure is a signal difference between light- and heavy-quarks: 
for light-quarks, Mf /irif > 20, while, for heavy-quarks, Mf /rrif < 5. This is 
a manifestation of the fact that the heavy-quark mass functions are well rep- 
resented by a constant whereas those of light-quarks are not; the accuracy of 
this approximation improving with increasing current-quark mass. The same 
is true of Zf{jP)\ i.e., Zj^lp"^) ^ 1. This feature provides the basis for a simpli- 
flcation of our study of semileptonic decays of heavy mesons: in the kinematic 
region explored by the decays one can approximate the dressed heavy-quark 



^ We employ a Euclidean space formulation with {'y^,^,^} = 25^^,, 7^ = 7^t and 
a ■ b = X)i=i C'ibi- A timelike vector, Q^, has < 0. 

^ Confinement entails that there is no real solution of + Mf{p'^) = 0; i.e., no 
"pole-mass" for the quarks. 



Sf{p) = -^7 ■ palip^) + (yiip^) = 
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Fig. 1. The behaviour of Mj(p^): u/d-quark - sohd hne; s-quark - long-dashed line; 
c-quark - dashed line; and 6-quark - dot-dash line, obtained using the model of 
Ref. [3]. See also Refs. [4,5]. The diagonal dotted line is p = M(p). 



propagator as 



S 



1 



f=c,b 



■ p + Mf 



(2) 



where Mf ^ M^, i.e., in the first instance, one can ignore the momentum de- 
pendence of gluon dressing of the heavy-quark propagator. In the DSE frame- 
work, "heavy-quark symmetry" [6] finds its foundation in this. 

Returning to the dressed quark- W-boson vertex, which describes the couphng 
of dressed-quarks to the W-boson, the vector piece satisfies the identity 

Q^^V^'^'{p; Q) = - V(P-) - {mf, - m;jrf ^^(p; Q) , (3) 



where = p + rj Q , p_ = p — i] Q , fj = (1 — rj) , and Q is the total momen- 
tum. Here rj is the momentum partitioning parameter that arises because of 
the arbitrariness in the definition of the relative momentum in a covariant 
formalism, and r{^-^^(p;Q) is the flavour-dependent scalar vertex. (In the ab- 
sence of interactions T{^'^'^{p;Q) = Ijj.) In studying the semileptonic decays 
of light mesons, the dressing of the vertex implied by this Ward-Takahashi 
identity was itself important, even neglecting the nonanalytic contributions 
mentioned in the introduction [1] . However, in the case of heavy-quarks, the 
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ability to neglect the momentum dependence of gluon dressing entails Eq. (2) 
and {rrif-^ — mf^)T^/^^ {p; Q) ^ (Mj^ — Mf^)Ij:,] and hence Eq. (3) is satisfied 
approximately by the bare vertex. This justifies the approximation 

K^^^^p; Q) = 7/. , (4) 



amplifying the heavy-quark simplification in the study of these decays. 
2. Semileptonic and leptonic decays. 

Herein we study the impulse approximation to the semileptonic D^i^u 
decay amplitude, defined by 



{D-{K)\h,c\B\P)) ^ + P), - f4t){K - P), (5) 

= N,j ^tro [r^- {^ + r)P; ~K) + r/(P + K))Tbo {I + f]K- P) x 
S^{t -7]P + f}K)iVl;^\l + (1 - 2'q)K- K - P)S^{£ + f}P- 7]K) 

where t = -{P - Kf, TBO,D-{k; Pf = C^TBO,D-{-k; P)C, C = 7274 is the 
charge conjugation matrix, and F^o ^j- are the meson Bethe-Salpeter ampli- 
tudesQ normalised via 



Nc I (0itrz) [Tbo,d- (k; -P)d';Sa{k + vP)Tb^.^d- (k; P)S„,{k - r/P) 

+ Tb^^d- {k; -P)Sd{k + fiP)TB^^D- {k; P)d^^S^/,{k - t^P)] = 2P^ , (6) 

which is the consistent, canonical normalisation in impulse approximation [7]. 

Simultaneously we study the leptonic decay of a heavy, pseudoscalar meson, 
M, which is described by one, dimensioned coupling, fm'- 



(0|/27M75/l|$Af(P))= (7) 

fMP, = N,j -0^^trn[l,l,SfAk + fiP)TM{k-,P)Sak-r^P)] . 

With this normalisation, /^r — 131 MeV. Equation (7) is exact but approxi- 
mation enters in the calculation of the Schwinger functions that appear in it; 
i.e., the quark propagators and Bethe-Salpeter amplitude. 



^ In solving a Bethe-Salpeter equation to obtain the bound state mass and F, the rj- 
dependence of F ensures that physical observables, such as the mass, are independent 

of T]. 
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2.1 Heavy-quark limit. At this point, no assumption has been made about the 
form of the Schwinger functions in Eqs. (5) and (7). In a study of hght-meson 
decays, these functions were determined from model DSEs, and the integrals 
then evaluated [1] . For heavy mesons, where the apphcation of DSEs has not 
hitherto been extensive, we will explore the consequences of Eqs. (2) and (4), 
in an estimation of the behaviour of the form factors and decay constants: an 
exploration of a heavy-quark limit in our framework. 

To be explicit, we set rj — 1 in Eqs. (5) and (7) (recall fj = 1 — rj), so that the 
heavy-quark carries all the momentum of the heavy meson, in which case the 
Bethe-Salpeter amplitude constrains the momentum of the light-quark in the 
meson. We write the heavy- meson momentum as = [Mf^ + E)v^, where 

E = Mfj — MfQ , with Mh the mass of the heavy meson, and is a timelike 
unit- vector, v"^ = —1. Equation (2) then yields 




(8) 



The bound state amplitudes, present in the integrands that arise in the cal- 
culation of physical processes, ensure that |p|/M/q <C 1. 

Using Eq. (2), or (8), the dressed quark- W-boson vertex is given by Eq. (4). 
For the heavy-meson Bethe-Salpeter amplitude we employ the Ansatz 

rso,n-{k;P)-i^{i-kh-v)^^, (9) 



which is motivated by the Bethe-Salpeter equation studies of Ref. [8]. This 
function describes the momentum distribution of the light-quark, which we 
assume to be the same in both B and D mesons. Our estimate of the leptonic 
decay constants and form factors will be robust if the results are insensitive 
to the form of <f{k'^). 

Using Eqs. (4), (8) and (9), the normalisation condition, Eq.(6), can be written 







where z — u — 2Ea/u and H — B^,D , . . . , etc., with the light-quark flavour, 
/, chosen appropriately. Clearly J^f^fMj^f — constant = -^^/-^ij/- Similarly, 



5 



we find for a lieavy-meson leptonic decay constant: 



/h = ^^ / duiV^-E)^{z) {aiiz) + kV^4i^)) , (11) 



whicli entails tliat 



fny Mh = constant , (12) 



except for deviations due to differences between u/d- and s-quark propagation 
cliaracteristics. 

Considering Eq. (5), we find, at leading order in 1/Mh, 

1 1 °° 

= ''^32^ I'^'^W I {""sM + ^^M) , (14) 





witli W = 1 + 2r(l -t){w- 1), zw = u- 2EJu/W and □ 



w 



2MbMd 



vb-vd ■ 



(15) 



Tfie canonical normalisation of the Bethe-Salpeter amplitude automatically 
ensures that ^(w = 1) = 1. Equation (13) is an example of a general result 
that, in the "heavy-quark limit" the semileptonic decays of heavy mesons are 
described by a single, universal function [6]. 

Equation (14) yields the result: > Pmin; where0 



2 ^ d^^jw) 
^■""^ ~ dw 




o y 2/o~rfMVP(u)2 [aUu) + ^at.{u] 



(16) 



The minimum physical value of w is tOmin = !> which corresponds to maximum 
momentum transfer with the final state meson at rest; the maximum value is ifmax — 
(M^ + M1))/{2MbM£)) = 1.6, which corresponds to maximum recoil of the final 
state meson with the charged lepton at rest. 

^ For (p{z) and ay^g{z) non-negative, non-increasing, convex-up functions of their 
argument, which includes ip = constant and a free-particle propagator, —CBiw) takes 
its minimum value at E' = 0. 
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which entails an upper bound on the slope of C{w): ^ |, independent of 
the details of our model and consistent with Ref. [9]. 

2.1 Light-quark propagator. The leptonic decay constant and ^(w) depend 
on the light-quark propagator and the momentum distribution of the light- 
quark in the heavy meson, described by fik"^)- Many observables involving 
hadrons composed of light-quarks have been studied using the u/d- and s- 
quark propagators specified by 



4{x) = 2%jr(2(x + fh})) + J'{hix)T{hx) {hi + hiT{ex)) (17) 



where ^{ij) = {I - c-y)/y, x = p'^/{2D) and: a{,{x) = 2Da(.{p^); a{.{x) = 
y/2D (Jgip^)] rhf = mf/y/2D, with D a mass scale. This form is motivated 
by extensive studies of the DSE for the dressed-quark propagator [2] and 
combines the effects of confinement and dynamical chiral symmetry breaking 
with free-particle behaviour at large spacelike-p^; i.e., asymptotic freedom. 

The parameters m^, bL, in Eqs. (17), (18) were determined in a x^-ht to a 
range of hght-hadron observables, which is described in Ref. [10] and leads to 
the values in Table 1. In the fit, the difference between the u- and d-quarks was 
neglected and only that minimal difference between u- and s-quarks allowed 
that was necessary to ensure: (ss) < (uu); and rrLs/mu ^ 1. (e = 10""^ is 
included in Eqs. (17), (18) only for the purpose of separating the small- and 
intermediate-p^ behaviour of the algebraic form, characterised by bo and 62; 
a separation in magnitude observed in numerical studies.) These light-quark 
propagators have since been used to successfully explore and predict a range 
of other hadronic observables [1,11]. We use them herein. 

Table 1 

The values of bf 3 are underlined to indicate that the constraints 6f 3 = 6" 3 were 
imposed in the fitting. The scale parameter D = 0.160 GeV^. 



m/ 


b^ 


b( 


b^ 


b^ 


u 0.00897 


0.131 


2.90 


0.603 


0.185 


s 0.224 


0.105 


2.90 


0.740 


0.185 



3. Calculated heavy- meson observables. 

In this study we have introduced one free parameter, E = Mh — M^jj^, which 
is a measure of the binding energy in a heavy meson, and an Ansatz for the 
Bethe-Salpeter amplitude, Eq. (9), which involves a single function, <^(/c^). 
Our calculated form of C{w), and the heavy-meson leptonic decay constants. 
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will depend on these quantities. Our predictions will be robust if they are 
independent of the details of our Ansatz for (f{k'^). 

To explore this we considered the following four one-parameter forms: 

(PAik"^) (Psik'^) fc{k'^) (Poik'^) 



exp - 



A^ 



A^ 



A2/ A;2 + A2 \k^ + A^ 



e 1- 



k^ 

A2 



(19) 



For a particular form, to fit the parameters (-E, A), we performed a fit to 
the following parametrisation of the experimental data [12] 



w + 



- ex.p ( (1 - 2p'^)- — -] , p2 = 1.53 ±0.36 ±0.14, 
1 V w + 1J 



(20) 



to fo = 216 ± 15MeV and /s = 206 ± 30MeV, which, in the absence of 
experimental data, is our weighted average of the lattice-QCD results collected 
in Ref. [13]. Using Md = 1.87 GeV, Md, = 1.97 GeV and Mb = 5.27 GcV, we 
obtain the results presented in Fig. 2 and Eq. (21), energies in GeV and 
dimensionless. 





E A 


Id fos fs 


A 


0.640 1.03 


0.227 0.245 0.135 1.55 


B 


0.567 0.843 


0.227 0.239 0.135 1.56 


C 


0.612 1.32 


0.227 0.242 0.135 1.55 


D 


0.643 1.02 


0.272 0.296 0.162 1.21 



(21) 



4. Observations and Conclusions. 

One observes from Eq. (21) that, independent of v^(/c^), a description of the 
data requires a "binding energy" in B, D mesons of ~ 0.6 GeV, which is 
consistent with the values obtained in Bethe-Salpeter equation studies of these 
systems: £^ ~ 0.4 — 0.6 [5,14]. Further, the size-parameter A ~ 1 GeV; and 
comparing this with an analogous parametrisation of the pion Bethe-Salpeter 
amplitude shows A ~ (1.5-2.0) A^r [15]. Hence, our Bethe-Salpeter amplitude 
represents the heavy meson as on object of small spacetime extent; occupying 
an intrinsic volume of only 5-20% that of the pion. 

Globally, considering Eq. (21) and Fig. 2, the Ansatze: A,B,C, for the Bethe- 
Salpeter amplitude provide equivalent descriptions of all observables. The step- 
function, Ansatz D, which is qualitatively different in form and corresponds 
simply to the introduction of a cut-off in the integrals that yield observables. 
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w 

Fig. 2. Our calculated form of S,{w) compared with experimental fits [12]. Experi- 
ment: dashed line, Eq. (20); short-dashed line, the linear fit ^{w) = \ — p'^{w — 1), 
= 0.91 zb 0.15 zb 0.06. Calculation: solid line, ^{w) calculated with forms A, B or 
C in Eq. (19) using the parameters in Eq. (21), there is no discernible difference on 
the scale of this figure; dot-dash line, ^{w) calculated with form D. 

not surprisingly leads to a poorer description. The fact that even this form is 
not completely unreasonable is due to the presence of the light-quark prop- 
agators in our study, which describe the confinement of light-quarks in the 
heavy meson and restrict their propagation range. It is the combination of the 
light-quark propagators, determined in studies of light-hadron observables, 
and the heavy-meson Bethe-Salpeter amplitudes that prescribe the values of 
the observables; i.e., as observed in Ref. [16], these observables are strongly 
influenced by the nonperturbative characteristics of confinement, dynamical 
chiral symmetry breaking and bound state formation. 

In our simple analysis, the calculated values of the leptonic decay constants 
are in quantitative agreement with the results in Ref. [5]. In particular, with 
the reahstic Ansatze: A,B,C, we note that fosl Id ~ 1-07, whereas Eq. (12) 
suggests the result fojlo = ^Md/Md, = 0.97. This difference, ~ 10%, is 
an indication of the influence and magnitude of light-quark effects in heavy- 
meson observables. Further, we note that while our value of fo = 227 MeV, 
agrees with that found in numerical simulations of lattice-QCD, our value 
of /b — 135 MeV is ~ 35% lower, and there appears no way to reduce this 
discrepancy. Our quantitative agreement with Ref. [5] suggests that the value 
inferred from lattice simulations may be erroneous. 

We obtain significant curvature in ^(w) on the physically accessible domain 
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in B — i> D^v decays: 1 < w < 1.6; i.e., as illustrated in Fig. 2, a comparison 
of our calculated result with the linear fit of Ref. [12] is unfavourable. In 
the present study, no choice of our two parameters yields a straight line, and 
trying to fit the short-dashed line in Fig. 2 leads to a value of ^ (wmax) ~ 20% 
higher than that inferred from experiment. We note that using any of the 
three, realistic, one-parameter Ansatze for the Bethe-Salpeter amplitude, our 
calculated result is in exact agreement with the experimental fit to the data, 
Eq. (20). The agreement between our calculation and the other nonlinear 
fits described in Ref. [12]: exp(-p2(w - 1)), = 1.27 ± 0.29 ± 0.12; and 
{2/{w + l)Yf\ p2 = 1.42±0.32±0.13, is not quite as good, with the deviation 
reaching ~ 5% at Wmax- From this we conclude that in the range ~ 1.2- 1.6 
is admissible, with ^ 1.55 most likely. 

Our calculation is not alone in indicating significant curvature in ^(w); there 
are qualitative similarities with the results of Refs. [16-18], for example. How- 
ever, our calculation is unique in being able to encompass a value of ^(wmax) ^ 
0.5 with physically reasonable and internally consistent values of our two pa- 
rameters, E and A. The detailed, constrained description of the light-quark 
component of the heavy meson, via its propagator and its momentum distri- 
bution function (Bethe-Salpeter amplitude), is responsible for this. 

A quantitative step in extending this study is a thorough analysis of the DSE 
for a heavy-quark, such as being explored in Ref. [19], which could provide the 
basis for an estimation of corrections to the "heavy-quark" limit expressed by 
Eqs. (2) and (8). One could then employ the calculated heavy-quark propaga- 
tor in Bethe-Salpeter equation studies of heavy-meson bound states, perhaps 
in a manner analogous to Ref. [20]. In this way E and A, which are fitting- 
parameters herein, would be correlated and their values fixed by the mass scale 
in the gluon propagator. This would lead to predictions of leptonic decay con- 
stants and semileptonic form factors that are less dependent on existing data. 
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